Inclusive 7r and K$ in 77 reactions at L3 
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Abstract. The e + e~ — ► e + e _ 7r° + X and e + e~ — > e + e~Kg + X reactions are studied 
at LEP using data collected with the L3 detector at y/s — 189 — 202 GeV. Preliminary 
results for the differential cross section da/dpx in the transverse momentum range 
fN| \ 0.2 GeV < px < 7.5 GeV at central values of pseudo-rapidity, \r]\ < 0.5 or | | < 1.5, 

and the differential cross section da/d\rj\ for different values of px are presented. For 
Pt < 1 GeV, the pt— dependence of the data are well described by an exponential fit. 
For px > 1 GeV the cross sections are compared to NLO QCD calculations and to LO 
^-j. \ Monte Carlo predictions. 
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INTRODUCTION 

Two-photon collisions are the main source of hadrons at LEP via the process 
e + e~ — > e + e~7*7* — > e + e~ hadrons . In the Vector Dominance Model (VDM), each 
photon can fluctuate into a vector meson with the same quantum numbers as the 
photon, thus initiating a strong interaction process with characteristics similar to 
^ . hadron— hadron interactions. This process dominates in the "soft" region, where 

hadrons are produced with a low transverse momentum p-r . Hadrons with high 
Pt are produced by the QED direct process, 7*7* — > qq, or by QCD processes 
originating from the hadronic content of the photon. Depending on wether one 



or two photons fluctuate into a hadronic system, the QCD processes are called 
single— resolved or double— resolved respectively. In two-photon interactions, QCD 
calculations are available for single particle inclusive production at Next-to-leading 
order (NLO) [1]. 

In this paper the inclusive 7r° production and the inclusive K s production from 
quasi— real photons are measured. The center-of-mass energy of the two interacting 
photons, W^, is greater than 5 GeV for 7r° and greater than 10 GeV for K s . The 
7T° is measured in the transverse momentum range 0.2 < px < 7.5 GeV and in the 
pseudo— rapidity interval —.5 < i] < .5, with 77 = — ln(tan#/2). The angle 6 is the 
polar angle of the 7r° , relative to the beam axis. The K s is measured in the range 
0.2 < p T < 4 GeV and -1.5 < 77 < 1.5 

The differential cross sections are compared to the Monte Carlo models PHOJET 
[2] and PYTHIA [3] and to analytical NLO QCD calculations [1]. 



EVENT SELECTION 



The data used for this analysis were collected by the L3 detector [4] in 1998 and 
1999 at beam energies from 189 GeV to 202 GeV, with an average energy value of 
194 GeV. The integrated luminosity is 414 pb _1 . 

The e + e~ — > e + e~ hadrons processes are simulated with the PHOJET 1 and 
PYTHIA 2 event generators. For the background the annihilation processes e + e~ — > 
hadrons^), ZZ( 7 ), Zee( 7 ), Wez/( 7 ) are simulated with PYTHIA [3]; KORALZ [5] 
is used for e+e" -> r+r-( 7 ) and KORALW [6] for e+e" -> W+W". For the 
e + e~ — > e + e~r + T~ channel the generator DIAG36 [7] is used. The events are 
simulated in the L3 detector using the GEANT [8] and GEISHA [9] programs and 
passed through the same reconstruction program as the data. 

The selection of e + e~ — > e + e~ hadrons events is based on information from 
the central tracking detectors and from the electromagnetic (BGO) and hadronic 
calorimeters. The following cuts are applied: 

1. To reject annihilation events, the total energy in the calorimeters is required 
to be less than 40 % of the center-of-mass energy. 

2. To exclude radiative events of the type e + e~ — > 7Z — > 7 qq, the total energy 
in the electromagnetic calorimeter must be lower than 50 GeV. 

3. Hadronic events are selected with at least 6 particles in the detector. These 
particles can be charged tracks or isolated clusters in the calorimeters. 

4. Quasi— real two-photon interactions are selected by excluding events with a 
cluster of energy greater than 70 GeV and polar angle greater than 33 mrad. 

5. Beam-gas and beam- wall events are suppressed by the requirement that at 
least 500 MeV must be deposited in the electromagnetic calorimeter. 

6. The analysis is limited to events with a visible hadronic mass W vis > 5 GeV , to 
allow comparisons with the Monte Carlo generators which generate events with 
Wjj> 3 GeV. Wy\s, is the effective mass of the hadronic system, calculated 
from the four- momentum vector of all the measured particles 3 . 

After these cuts, the background (~ 1 %) is due to e + e~ — > hadrons and e + e _ — > 
e + e~r + T~ events. 

Inclusive 7r° production is studied via the n° decay into two resolved photons. A 
photon is defined as an electromagnetic cluster formed by the energy deposited in 
at least 2 BGO crystals, with energy greater than 100 MeV and separated by more 
than 200 mrad from the closest track. The number of photons detected is required 
to be at least 2 but not more than 20. 

V PHOJET versionl. 05c 

2) PYTHIA version 5.718 and JETSET version 7.408 

3 ) All particles are considered to be pions, except for isolated clusters in the electromagnetic 
calorimeter which are considered to be photons. 



Almost 2 10 6 events are selected, leading to more than 8 10 7 two photon com- 
binations. The distribution of the mass of the reconstructed 77 system shows a 
narrow 7r° peak. A fit of this peak with a gaussian function over a Chebyshev 
polynomial background gives a resolution of 7.29 ±0.03 MeV. 

Inclusive Kg production is studied via the Kg decay into n + n~ . The distance, 
in the transverse plane, between the secondary vertex and the e + e~ interaction 
point is required to be greater than 3 mm. The angle between the flight direction 
of the K s candidate (taken as that of the line between the interaction point and 
the secondary vertex in the transverse plane) and the total transverse momentum 
vector of the two outgoing tracks must be less than 0.075 rad. After this cut, 
around 5 10 5 events are selected. A fit of the K s peak gives a resolution of 10 ± 1 
MeV. 



DATA ANALYSIS 

To evaluate the number of 7r°s or K®s, on each bin of and \rj\, a fit is performed 
using a gaussian for the signal and a Chebyshev polynomial parametrisation of the 
background. In the n° case, to estimate the uncertainty on the background subtrac- 
tion a side-band background subtraction is also performed. A three sigma region, 
"the signal region", is defined around the 7r°peak. Two "background regions" are 
defined on each side of the signal region with the same width. The average number 
of 77 combinations in the background regions is subtracted from the number of 
combinations in the signal region. Both methods give consistent results. 

The reconstruction efficiency is evaluated using the two Monte-Carlo generators: 
PHOJET and PYTHIA. The kinematic limits of the cross section measurement are 
defined by imposing, at the generator level, > 5 GeV for 7r° and W 77 > 10 GeV 
for Kg . The virtuality of the photon is limited to Q 2 < 1 GeV 2 in PYTHIA, in 
PHOJET, a limit Q 2 < 3 GeV 2 for tt° and Q 2 < 8 GeV 2 for K° s is imposed. We 
have verified with PHOJET that the results do not depend on the Q 2 cutoff. The 
number of it / K s reconstructed is obtained with the same procedure used for the 
data. As the two generators reproduce equally well the experimental distributions 
[10], a weighted average of the two is used to correct the data. 

The level— 1 trigger efficiency is evaluated by comparing the response of two inde- 
pendent triggers: the energy trigger [11] and the track trigger [12]. The efficiency is 
90—95 %, varying with the data taking conditions during the year. The efficiencies 
of higher level triggers, calculated using prescaled events, vary from 85% at low p T 
to 100% at high p T . 

The main uncertainty on the measured cross sections comes from the choice of 
Monte Carlo generator used to calculate the reconstruction efficiency and from 
the background estimation. Half of the difference between the two generators and 
between the two estimations of the background are taken as systematic errors. The 
uncertainties due to cut variations are negligible. 



Statistical and systematic errors are added in quadrature in the following. All 
the results are preliminary. 

RESULTS 

The differential cross sections da / dpx are shown in fig. 1. Exponential fits 
e - PT /<PT> f or 0.2 GeV < p T < 1.0 GeV reproduce well the data with < p T >= 
233 ± 1 MeV for ir° and <pr >= 357 ± 9 MeV for Kg . This behaviour is charac- 
teristic of hadrons produced by soft interactions and the fitted values of <pr> are 
similar to the one obtained in hadron— hadron or photon— hadron collisions [13]. 

Due to the existence of the direct process and of hard QCD interactions, 7*7* col- 
lisions start to exhibit a higher cross section at px values greater than ~ 1.5 GeV. 
In fig. 2 and 3 the da/dpT&mi da/d \rj\ differential cross sections are compared to 
analytical NLO QCD predictions [1]. For this calculation, the flux of quasi— real 
photons is obtained using the Equivalent Photon Approximation, taking into ac- 
count both transverse and longitudinal virtual photons. The interacting particles 
can be photons or partons from the 7 —>■ qq quantum fluctuation, which evolves, 
via the Altarelli-Parisi equation, into quarks and gluons. The Gordon and Storrow 
[14] parton density functions are used. All elementary 2^2 and 2^3 processes 
are considered. The renormalization scale, the factorisation scale and the fragmen- 
tation scale are taken to be equal: \i — M — Mp = ^p?. The uncertainty in the 
NLO calculation is estimated by varying the value of £ from 0.5 to 2.0. 

The measured cross sections are one order of magnitude above the direct process 
predictions. In the ir° case, for px > 2 GeV, the data are also higher than the 
predictions of NLO QCD calculations. 
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FIGURE 1. The differential cross sections da/dpr as a function of pr for ir° and Kg . The low 
Pt spectrum is well reproduced by an exponential fit, characteristic of soft interactions. 
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FIGURE 2. Comparison of the da/dpr distributions with NLO-QCD calculations. The central 
values (full lines) are the calculations with £ = 1 (see text), the dashed lines refer to £ = 0.5 (upper 
lines) and to £ = 2 (lower lines). The contribution of the direct process alone are indicated with 
dash-dotted lines. 



The cross sections are also compared to Monte Carlo predictions in fig. 4. The 
high pt region is better reproduced by PYTHIA than by PHO JET in the 7r° case. 
For the K s both Monte-Carlos reproduce quite well the data. 

The differential cross sections, calculated for W 77 > 10 GeV, are compared to 
OPAL results [15] in fig. 5. OPAL measured the inclusive production of charged 
hadrons (mainly tt^) and of Kg in the range \rj\ < 1.5. The experiments agree 
within the experimental errors. 




FIGURE 3. Comparison of the da /d\ -q | distributions with NLO-QCD calculations. The lines 
are as in fig. 2. 
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FIGURE 4. Comparison of the da/dpT distributions with Monte Carlo predictions. 
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FIGURE 5. The L3 and OPAL da/dpr distributions compared. The OPAL data are scaled to 
the same center-of-mass energy and, in the n° vs. charged hadron case, are scaled also to the 
same pseudorapidity range and fragmentation function. 



CONCLUSIONS 



This analysis is the first study of inclusive 7r° production in two photon collisions 
at LEP. At low p T , the differential cross sections da /dpx of ir° and K s production 
follow an exponential law similar to that previous observed in hadron-hadron and 
hadron-photon collisions. At higher px , the contributions of 77 — > qq and hard 
QCD processes are clearly visible. The agreement with previous measurement is 
very good, but there is only a qualitative agreement with the predictions of LO 
Monte Carlo and NLO QCD calculations. 
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